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PREFACE 


\ 


The  development,  implementation,  test  and  evaluation  of  an 
interactive  FORTRAN  program  for  Covariance  analysis  was  part  of 
the  GPS  software  development  project  at  DMAHTC,  Geodesy  and 
Surveys  Department,  Techniques  Office  (GST). 

It  .was  performed  under  the  supervision  of  Ben  Roth,  Chief, 
GST  and  Professor  Dhaneshwar  Hajela,  Department  of  Geodetic 
Science,  The  Ohio  State  University. 

This  report  was  submitted  to  the  Graduate  School  of  The  Ohio 
State  University  as  partial  fulfillment  of  the  requirements  for 
the  Master  of  Science  degree  in  Geodetic  Science. 
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ABSTRACT 


A  transportable,  generic,  interactive  FORTRAN  program  for  the 
covariance  analysis  of  pseudo-range  experiments  in  the  Global 
Positioning  System  (GPS)  is  presented  for  use  in  simulation 
studies,  experiment  planning  and  systems  development.  The 
interactive  mode  is  especially  suited  for  covariance  analyses 
involving  multiple  parameters,  by  allowing  ease  of  operation  for 
data  production  and  providing  savings  in  time  when  computing 
geodetic  point  positioning. 

The  mathematical  models  (presented  in  Chapter  3)  on  which  the 
program  is  based  provide  an  overview  of  the  GPS  point  pos¬ 
itioning  process,  emphasising  the  areas  of  interest  to  geodesy. 
Also,  special  emphasis  is  placed  upon  station  tracking  geometry. 
The  FORTRAN  program  allows  for  understanding  correlation 
coefficients  between  the  various  parameters.  A  separate  model 
suitable  for  covariance  analysis  is  also  presented. 

This  interactive  FORTRAN  program  has  been  implemented,  on  the 
UNIVAC  at  DMAHTC.  The  program  was  tested  with  simulated  GPS 
data.  The  results  indicate  that  the  program  is  effective  and 
efficient  in  helping  to  establish  the  reliability  of  the  GPS 
point  positioning  technique.  The  program  can  be  used  in  Simula-* 
tion  studies,  planning  test  and  evaluation  experiments,  and  GPS 
software  development.  The  quest ion/ answer  method  has  been  used, 
to  develop  this  interactive  program  which  enables  the  user  to 
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1.  INTRODUCTION 


"...Many  users  of  satellite  point  positioning  have 
adopted  a  wait-and-see  attitude  toward  GPS..." 

(Brown,  1979). 


1.1  Background 

The  Defense  Mapping  Agency  Hydrographic/Topographic 
Center  (DMAHTC)  provides  mapping,  charting  and  geodetic  (MC&G) 
support  to  the  Department  of  Defense  (DoD)  and  other  federal 
agencies. 

The  Geodesy  and  Surveys  Department  (GSD),  DMAHTC,  has  thm 
mission  to  acquire  and  reduce  satellite  tracking  data  to  support 
MC&G  requirements. 

Currently,  GSD  utilizes  the  Navy  Navigation  Satellite  System 
(NMSS)  to  support  MC&G  requirements.  This  utilization  includes* 
operation  and  maintenance  of  a  worldwide  satellite  tracking 
network  (TRANET);  data  communication,  processing  and  computation 
of  the  precise  ephemerides  and  the  deployment,  operation  and 
maintenance  of  portable  satellite  receivers  (Geoceivers  and  MX 
1502s)  for  the  determination  of  precise  geodetic  positioning 
points  and  relative  positions. 


In  the  future,  6SD  will  switch  to  the  NAVSTAR  Global 
Positioning  System  to  support  MC&G  requirements  when  it  replaces 
NNSS  in  the  late  1990s  as  OoD's  prime  satellite  navigation 


system. 


System  Description 


A  description  of  the  system  should  begin  with  the  NAV5TAR 
Global  Positioning  System  (GPS).  GPS  is  a  space  based  radio 
navigation  system  under  development  by  DoD,  GPS  Joint  Program 
Office  (JPO).  The  system  is  being  developed  to  support  the 
military  forces'  need  to  precisely  compute  their  position, 
velocity  and  time  any  place  on  or  near  the  earth,  24  hours  of 
every  day,  under  all  weather  conditions,  and  within  a  common 
reference  frame. 

GPS  completed  the  concept  validation  (Phase  I)  in  1979 
and  is  currently  in  full-scale  development  (Phase  II).  The 
system  development  schedule  calls  for  an  initial  operational 
capability  (IOC)  in  the  1986-87  time  frame  and  a  full  opera¬ 
tional  capability  (FOC)  in  the  1988-1989  time  frame.  At  FOC,  GPS 
will  replace  NNSS  as  DoD's  prime  satellite  navigation  system. 

GPS  consists  of  3  parts,  a  Space  Segment,  a  Control 
Segment  and  a  User  Segment.  Dfard,  1981]. 

Space  Segment 


The  operational  Space  Segment  will  consist  of  18  satel¬ 
lites  configured  into  6  equally  spaced  55  degree  inclination 


orbit  planes.  Each  plane  will  contain  3  equally  spaced  satel¬ 
lites  which  are  each  in  12  hour,  nearly  circular  (26,000 
kilometer  radius)  orbits.  The  resulting  satellite  constellation 
will  provide  on  or  near  the  earth,  continuous  user  visibility 
to  a  minimum  of  4  satellites.  In  addition  to  the  18  satellites, 
a  number  (to  be  determined)  of  spare  satellites  will  be  in 
orbit  to  support  continued  FOC  in  the  event  of  failures.  The 
launch  of  the  operational  satellites  will  be  accomplished  with 
the  Space  Shuttle  commencing  in  1986. 

The  Space  Segment  of  Phase  I  and  Phase  II  comprises  4-5 
developmental  satellites  which  are  in  12  hour,  nearly  circular 
(26,000  kilometer  radius),  63  degree  inclination  orbits.  The 
satellites  are  configured  to  provide  4  satellite  visibility  at 
Yuma  Proving  Ground,  Arizona  (YPG) ,  which  is  the  principal  test 
location  for  the  User  Segment. 

The  GPS  satellite  contains  a  S  band  receiver,  a  L  band 
transmitter,  a  solar/battery  power  system,  a  signal  proces¬ 
sor/generator  and  a  redundant  system  of  precise  frequency/ time 
standards.  The  satellite  continuously  broadcasts  spread  spectrum 
signals  at  center  line  frequencies  of  1575.42  MHz  (LI)  and 
1227.60  MHz  (L2).  These  signals  are  encoded  with  the  satellite's 
ephemerides  and  clock  estimates,  which  are  transmitted  daily  to 
the  satellite  from  the  Control  Segment. 

Control  Segment 


The  Operational  Control  Segment  will  comprise  five 
monitor  stations  (MSs),  three  upload  stations  and  an  operational 
control  station  (OCS).  The  MSs  will  be  located  worldwide  and 


will  track  the  GPS  satellites  obtaining  measurements  of  range 
and  integrated  Doppler  at  LI  and  L2.  These  data  will  be  communi¬ 
cated  to  the  OCS  where  they  will  be  processed  to  provide  precise 
predictions  of  the  satellites'  ephemerides  and  clock  estimates. 
These  predictions  will  then  be  communicated  to  the  upload 
stations  where  they  will  be  transmitted  to  the  satellites  via 
S-band  for  broadcast  to  the  User  Segment. 

The  Control  Segment  of  Phase  I  and  Phase  II  is  comprised 
of  four  MSs,  one  upload  station  and  a  Master  Control  Station 
( MCS ) .  The  upload  station,  MCS  and  an  MS  are  collocated  at 
Vandenberg  AFB,  California.  The  remaining  three  MSs  are  located 
in  Alaska,  Hawaii  and  Guam.  This  configuration  supports  User 
Segment  testing  at  YPG. 

The  development  and  implementation  of  the  Operational 
Control  Segment  will  be  phased.  The  operational  MSs  will  be 
deployed  in  late  1984  and  interfaced  with  an  Intermediate 
Control  Station  (ICS)  which  will  be  collocated  with  the  develop¬ 
mental  MCS  at  Vandenberg  AFB.  In  1986,  the  OCS  will  begin  to  be 
phased  in  and  will  be  located  at  a  site  in  Colorado  along  with 
an  upload  and  monitor  station.  In  1987,  the  OCS  will  take  over 
the  system  and  ICS  will  remain  as  a  backup  system. 

User  Segment 

The  operational  User  Segment  consists  of  a  family  of  user 
equipment  (US)  which  will  support  a  variety  of  service-peculiar 
navigation  requirements.  The  user  equipment  will  comprise 
antennae,  receivers,  navigation  processors  and  flexible  modular 
interfaces  (FMIs).  The  FMIs  support  the  integration  of  the  UE 
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The  major  objective  of  this  study  is  to  present  a 
variance-covariance  analysis  for  use  in  systems  development, 
simulating  and  planning  geodetic  point  positioning  using  NAVSTAR 
Global  Positioning  System  (GPS).  An  interactive  program  is 
offered  with  an  explanation  of  the  theory  and  mathematical 
models  on  which  the  program  is  based.  An  overview  of  GPS  for 
those  interested  in  applying  the  pseudo-range  method  for 
geodetic  activities  is  also  presented. 

The  variance-covariance  study  is  based,  on  the  method  of 
estimate-consider  analysis  of  parameters.  The  "best"  estimation 
of  a  vector  of  unknown  parameters,  X,  usually  computed  via  a 
formulation  of  the  state  of  the  system,  called  estimate  parame¬ 
ters,  these  being  the  latitude,  longitude  and  height  of  a 
station  in  a  common  reference  f reuse.  Thus,  one  must  "consider" 
the  effect  of  the  unestimated  parameters,  called  consider 
parameters,  which  are  tropospheric  correction,  the  satellite's 
long-track,  cross-track  and  radial  track,  and  the  bias,  drift, 
and  aging  of  the  combined  user-satellite  clock  parameters.  The 
program  is  written  in  the  interactive  mode  to  allow  the  user  to 
input  the  necessary  estimate  and  consider  parameters  for 
analysis  of  the  standard  deviations  of  each  plus  the  correlation 
factor  between  each  parameter.  The  consider  covariance  shows  the 
effect  of  the  consider  parameter  on  the  estimate  parameter 


[Bierman,  1977].  Plots  are  supplied  with  the  simulation  study 
indicating  optimal  satellite  tracking  time  for  the  establish— 
ment  of  the  geodetic  point  using  GPS  point-positioning.  Also, 
the  program  allows  the  user  to  plot  the  standard  deviation  of 
estimate  parameters  and  correlation  coefficients  between  each 
estimate  parameter. 

The  choice  of  the  parameter  set  was  influenced  by 
previous  studies  by  Milliken  and  Zoller  (1980)  and  Fell  (1980). 
The  primary  emphasis  is  placed  on  the  effect  of  consider 
parameters  on  the  estimate  parameters  of  a  geodetic  point,  (see 
list  of  acronyms) 

The  introduction  of  an  interactive  program  for  the 
variance-covariance  study  allows  the  user  for  the  first  time  to 
simulate  a  geodetic  positioning  experiment  and  view  the  results 
in  real  time.  It  was  decided  to  write  the  program  in  a  generic 
sense  for  future  use  since  the  mathematical  models  can  easily  be 
expanded  or  changed. 

This  study  also  illustrates  a  technique  which  allows  the 
geodesist  to  analyze'  station  selection  for  optimal  network 
development.  It  provides  the  geodesist  with  the  tools  for 
understanding  the  process  of  measurement  and  the  capability  to 
partake  in  future  simulations  of  experiments  and  planning  for 
geodetic  point  positioning.  With  this  information  the  geodesist 
can  address  the  problem  areas  and  implement  sound  statistical 
analysis  of  data  and  development  of  improved  mathematical 
models. 
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Organization  and  Scope  of  Stud 


Chapter  2  covers  the  basic  geometry  of  ground  station 
tracking  and  a  study  of  Geometric  Dilution  of  Precision  (GDOP). 
The  description  of  this  geometric  study  involves  the  mathematics 
for  calculation  of  visibility  circles,  GPS  sub-satellite  points. 
Position  Dilution  of  Precision  (PDOP)  and  optimal  ground  station 
tracking  geometry  for  any  station  in  the  world.  The  value  of 
GDOP  itself  is  a  composite  measure  that  reflects  the  influence 
of  satellite  geometry  on  the  combined  accuracy  of  the  estimate 
of  user  time  and  user  position.  The  value  of  PDOP  is  that  value 
that  reflects  the  influence  of  satellite  geometry  precluding  the 
estimate  of  user  time.  Graphs  are  supplied  showing  the  lowest 
GDOP  value  for  tracking  GPS  satellites  from  a  ground  station 
over  a  24  hour  period  indicating  optimal  ground  station  trackings 
for  a  GEOSTAR  GPS  user.  The  region  around  a  ground  station  where 
a  GPS  satellite  is  visible,  determined  by  station  coordinates 
and  the  semi-major  axis  of  satellite,  is  a  visibility  circle.  An 
explanation  of  the  interactive  program  for  future  use  in  studies 
involving  the  optimization  of  ground-station  tracking  of  GPS 
satellites  is  given. 

In  Chapter  3  the  mathematical  models  used  in  the  vari¬ 
ance-covariance  study  are  described  as  well  as  future  possible 
model  refinements.  A  model,  suitable  for  covariance  analyses,  is 
presented  for  determining  the  correlation  coefficients  between 
the  estimate  and  consider  parameters.  A  sensitivity  study  using 
estimate  and  consider  parameters  from  which  geodetic  parameters 


are  estimated  is  also  described.  A  summary  of  the  sensitivity 
study  is  discussed  and  illustrated  by  using  graphics  developed 
on  the  Hewlett  Packard  9845  Desktop  Computer. 


It  was  found  that  the  models  discussed  in  Chapter  3  do 
provide  solutions  to  geodetic  positioning  using  GPS  pseudo-range 
measurements.  These  results  are  discussed  in  Chapter  4  along 
with  recommendations  for  further  investigations.  The  last 
chapter,  chapter  5,  also  addresses  some  of  the  problems  involved 
in  the  development  and  use  of  the  models  which  will  be  created 
by  users  of  the  program. 

The  papercontains  the  following  graphs:  PDOP  and  GDOP, 
sub-satellite  point*  on  a  world  map,  visibility  circles  on  a 
world  map,  elevation  angle  vs  time,  azimuth  vs  elevation  angle, 
and  standard  deviations  of  parameters  and  correlation  coeffi¬ 
cients  among  parameters  vs  time.  Appendix  A  contains  a  listing 
of  variance-covariance  interactive  FORTRAN  program.  Appendix  B 
contains  a  sample  run  as  viewed  at  the  interactive  computer 
terminal  screen. 


2.  GEOMETRIC  SIMULATION  STUDY 


2. 1  Introduction 

One  of  the  main  reasons  for  the  development  of  NAV- 
STAR/Global  Positioning  System  (GPS)  is  to  provide  GPS  users 
with  a  precise  navigational  fix  at  any  point  anywhere  in  the 
world.  This  chapter  describes  a  rapid  method  for  performing  the 
geometric  simulation  of  station  tracking  for  any  point  on  earth. 
In  the  future  when  the  NAVSTAR/GPS  is  fully  operational,  the 
simulation  method  described  in  this  section  will  not  only 
provide  a  rapid  nagivational  fix  to  any  observer  at  any  point  in 
the  world,  but  also  allow  the  user  to  verify  that  the  geodetic 
receiver  has  selected  the  most  appropriate  satellites  to  ensure 
proper  data  acquisition. 

The  procedure  for  simulating  a  GPS  satellite  ephemeris  is 
discussed  first.  The  graphic  section  in  the  appendix  illustrates 
world-wide  NAVSTAR/GPS  sub-satellite  ground  track  points,  and 
the  visibility  circles  of  the  possible  tracking  stations  the 
Defense  Mapping  Agency  (DMA)  will  deploy  for  the  testing  of  the 
GEOSTAR  geodetic  positioning  receiver.  Secondly,  the  graphics 
section  shows  strength  of  the  geometry  for  the  NAVSTAR/GPS  cons¬ 
tellations  for  these  respective  tracking  stations.  In  doing  this 
analysis,  we  describe  the  principle  of  Geometric  Dilution  of 
Precision  (GDOP)  and  demonstrate  the  results  using  computer 
graphics.  Finally,  there  is  a  brief  discussion  of  the  GPS 
receiver,  GEOSTAR,  that  will  be  utilized  by  DMA  in  the  near 
future  for  geodetic  point  positioning. 


2.2 


GPS  Keplerian  Orbit  Generation 

We  describe  the  orbit  generation  program  in  this  section. 


The  first  step  in  the  procedure  is  the  location  of  sub-satellite 
points,  and  the  production  of  visibility  circles.  The*  Kepler ian 
elements  of  each  satellite  from  NSWC  and  the  beginning  and  the 
ending  epoch  for  the  orbit  generation  are  read  into  a  subroutine 
that  converts  the  NSWC  orbital  elements  into  the  traditional 
Keplerian  orbit  elements.  The  value  for  the  rate  of  change  for 
a,  e,  and  i  are  neglibly  small.  Then,  the  necessary  input  such 
as  GM,  ae,  inverse  flattening,  and  angular  velocity  of  the 
earth's  rotation,  <ue  ,  (WGS  72)  plus  the  six  Keplerian  elements 
are  also  read  into  the  front  part  of  the  program. 

The  epoch  for  the  orbit  generation  observations  was  1983, 
Day  17,  b*1,  0®,  0.0  sec.  The  earth  centered,  earth-fixed 
position  and  velocity  vectors  are  computed  by  the  program  for 
each  GPS  satellite  at  5  minute  intervals.  CMueller,  19823.  The 
program  is  interactive  so  that  the  user  can  easily  change  the 
time  span  and  number  of  intervals /unit  of  time  of  the  orbit 
generation.  The  program  has  been  tested  by  a  comparison  with 
NSWC  ephemerides,  and  found  to  be  within  millimeter  precision 
for  several  satellites  at  the  epoch  1983,  Day  17,  0*1  0®  0.0  sec. 

The  program  also  has  the  capability  for  updating  Keple— 
rian  elements  so  that  the  user  could  reasonably  generate 
ephemerides  without  current  NSWC  elements.  This  method  is 
outlined  here  and  it  may  be  used  if  a  simulation  study  needs  to 
be  done  for  point-positioning  tracking  projects  in  the  future. 


First,  one  must  update  the  argument  of  perigee,  <*>  , right 
ascension  of  ascending  node,®  ,  and  the  mean  anomaly.  The  value 
for  the  rate  of  change  for  a,  e,  and  i  are  negligibly  small.  The 
rate  of  change  for  these  elements,  due  to  the  second  zonal  coef¬ 
ficient  C20  in  the  earth's  gravitational  potential,  is.  given  by 
[KAULA  1966,  page  39] 


a 


3  g  c20  *| 
4(l-«2)3/2a2 


(3cos3i 


1) 


(2.2  -1) 
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» 
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3  a  C2Q  a* 

2(l-«2)2  a2 


cos  (  i) 


i 


3  a  C20 
4(l-«2)2  a2 


(l-5cos^i) 


(2.2  -2) 


(2.2  -3) 


where  ae  is  earth's  mean  equatorial  radius,  a,  e,  i,  n  are 
respectively  the  semi-major  axis,  eccentricity,  inclination,  and 
mean  motion  of  the  satellite  orbit. 

During  the  program  run,  the  program  asks  if  the  elements  need  to 
l  be  updated  to  a  certain  day;  and  then  performs  the  update  in  a 

subroutine  by  multiplying  the  day  change  by  the  rate  of  change 
l  for  each  Keplerian  element.  One  interesting  finding  was  the  lack 

"  of  contribution  given  by  the  mean  anomaly  to  the  ephemeris 

i  caused  by  the  perturbations  of  the  satellite.  The  way  this  came 

C 

K 
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\ 
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about  involved  the  need  for  a  requested  station  tracking 
geometry  and  study  to  verify  the  rise  time  of  satellites  with 
another  agencies  updated  ephemeris.  This  simulation  created  the 
first  opportunity  to  do  the  geometry  analysis  with  cross 
checking  the  data  utilizing  ui,n  and  <u,n,M  orbit  parameter  update. 
Interestingly,  the  geometry  calculated  using  ,  parameter 

update  was  more  correct  and  obviously  quicker  to  compute  than 
with  all  three. 

2.2.1  Sub-Satellite  Points 

The  procedure  for  calculation  of  sub-satellite  points  is 
performed  by  using  a  method  developed  by  B.  R.  Bowring.  [1976 
Survey  Review  No.  181] 

The  program  procedure  for  calculating  the  sub-satellite  points 
is  outlined  as  follows: 

1)  Input  X, Y, Z  of  satellite;  also  inverse  flattening, 
semi-major,  a,  and  semi-minor ,b,  axes  of  reference  ellipsoid  and 
a  test  condition  for  iteration  of  ellipsoidal  latitude, 

2)  Solve  for  ZQ  and  XQ  where 
ZQ  *  ZA 

XQ  *  (XA  +•  YA )l/2 

3)  Solve  for  reduced  latitude, 

3  -  DAIAN2(b*ZQ,a*XQ) 

4)  Set  iteration  to  zero  and  intermediate  value,  ( “store" ) , 


to  zero 


5)  Compute  tangent  of  PHI 

can  4>  •  (ZQ  +  e'2*b*(dsin  8)3)  /  (XQ  -  a*e2*(dcos  3)3) 

6)  Test  condition  of  PHI  for  iterative  process. 

if  [DABS  (PHI  -  STORE) .LT.TESTC.  OR  . ITER.GT.10]  Go  to 

Step  9 

If  not  continue  with  Step  7. 

7)  Put  store  equal  to  PHI 

STORE  *  PHI 

8)  Compute  BETA. 

S  -  0AXAN2[ (1  -  e2)1^2  *  TAM  t  ] 

Go  To  Step  4 

9)  Compute  *  ,  Geodetic  Longitude,  \  «  TAN(YA.XA) 

This  program  produces  a  file  for  each  satellite's  ground  track 
points  so  that  the  points  can  be  plotted  on  a  world  map  using  a 
Hewlett  Packard  9845  desktop  computer.  The  plots  are  found  in 
the  following  sections. In  this  manner,  one  can  obtain  graphics 
illustrating  the  optimal  time  for  observing  the  present  GPS 
constellations  for  any  ground  station.  (See  sample  plots  in  this 
chapter).  These  graphs  show  the  relative  positions  of  each 
satellite  in  the  current  five  satellite  configuration.  The 
azimuth  versus  elevation  angle  versus  time  plots  add  confidence 
to  my  PDOP  and  GDOF  measurements  because  one  can  see  that  as  the 
satellites  spread  out  over  the  station,  the  value  of  GDOP  and 
PDOP  decreases. 


.2.2.2  Visibility  Circl— 

In  this  section#  the  geometry  and  mathematics  of  GPS  viaibil— 
[  ity  circles  for  specific  ground  stations  is  presented.  Given  the 
geocentric  coordinates#  and  Re  the  mean  radius  of  the  earth# 
the  coordinates  convert  in  the  following  way: 


(1)  ZSTN  *  *  sin  4 

(2)  HEQ  •  Ra  *  cos  <j> 

(3)  Determine  values  for  XSTN.and  YSTN  ,  Thus  we  have, 
(See  Figure  2.2-1) 

XSTH  -  HEQ  *  cos  X 
TSTB  ■»  HEQ  *  sin  A 
ZSTIff  *►  H.e  *■  sin. 


Figure  2.2-1  Geocentric  Station  Coordinate  System 


therefore , 

(4)  s  •  AacsuKCV^*9^900  +  a.)) 


(2.2-8) 


a  -  180°  -  (90°  +  at  *  S)  (2.2-9) 

51  *  ^  *  sla  U  (2.2-10) 

ii  •  *  cos  &  (2.2-11) 

where  M  is  the  distance  along  the  station  position  vector 


and  N  is  the  distance  perpendicular  to  the  station  position 
vector  that  locates  a  point  on  the  visibility  circle. 

An  assumption  is  made  that  Rs  is  the  mean  semi  major  axis 
for  any  GPS  satellite  and  is  equal  to  approximately  26500 
km.  Another  assumption  is  that  the  elevation  angle  cutoff 
for  each  satellite  is  15  degrees. 

The  visibility  circle  is  then  determined  by 


HU 


Figure  2.2-3  Overhead  View  of  the  Visibility  Circle 


Figure  2.2—3  Polar  71 as  of  the  71slblU.tr  Circle 

First,,  by  solving:  for  x,and  y,  in:  the  following-  wayt 
r  »  It  *►  cos  ol 

T  »  IT  *  slit  a.  (See  Figure  2*. 2-3) 

*•  •  »  -  H^*tos  IT 

where  a.  is  fros  0°-  to  360°  in  1°  steps. 

After  computing-  the  necessary  quantities, .  one  plots  the  X 
and:  t  coordinates  tor  alpha.,,  cl  (see  Figure  2-3)  frost  0°  to  360° 
in-.  L  degree  steps*  creating:  the  visibility  circle-  It  i»  easily 
seen:  free  this  ezasple  that  one  can  rotate  the  north-  pole; 
situation  by  using-  *  and  \  ,  with  RL  and.  R3  rotation,  matrices. 
[Mueller ,  196 9 X  for  plotting:  visibility  circles:  for  any  station, 
in:  the  world- 


Transit  satellite  visibility  circle  vs.  GPS  satellite  visibility  circle 
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Thus,  from  the  visibility  circle  plots ,  one  has  an  idea  of 
the  tremendous  advantage  of  ground  station  tracking  GPS  has  over 
the  present  transit  navigation  system.  The  observer  has 

increased  time  of  observation  and  as  many  as  11  satellites  in 
the  1990 's  depending  on  user  location  and  orbit  inclination  to 
observe.  CMilliken  and  Zoller,  1980].  (See  Figure  2.2-4) 

2. 3  Strength  of  Geometry 

The  geometry,  which  is  formed  by  the  observations  between 
the  tracking  station  and  each  visible  GPS  satellite  determines 
how  the  measurement  errors  are  propagated  to  the  estimated 

parameters  in  the  Least  Squares  adjustment.  The  strength  of  the 
geometry  to  minimize  the  magnitude  of  the  error  propagation  is 
measured  by  the  Geodetic  Dilution  of  Precision  CMilliken  and 
Zoller  1980].  This  section  describes  the  technique  for  a  GDOP  & 
Position  Dilution  of  Precision  (PDOP)  solution. 

The  geometry  of  the  GPS  const el latione,  usually  four  selected, 
satellites,  contributes  to  the  magnitude  of  the  user  accuracy 
and  also  to  the  potential  for  position  errors  in  the  GPS 

navigation  fix.  Any  GPS  user  has  to  have  a  special  interest  in 

this  fact  because  the  four  "best"  satellites  selected  by  the 
user  receivers  are  .those  with  the  lowest  GDOP.  CMilliken  and 
Zoller  1980]. 

The  covariance  matrix  of  the  error,  Xu,  in  the  estimate  of 
Xu  is  given  by:  (1) 

°u  *  «&„>  1guCov  5<V  -  (2.3-1) 

GDOP  is  calculated  by  setting  the  covariance  matrix 


L8 


Cov  (Ay  -p)  equal  to  the  identity  matrix. 

Thus,  the  remaining  portions  of  (1)  can  then  be  reduced  to  the 
following: 


-  cgJ  Ou)-i 


(2. 3.-2) 


Cov  (Au  S-P)  primarily  reflects  range  measurement  error  statis¬ 
tics  (i.e.  satellite  ephemeris,  ionospheric  model,  and  instru¬ 
mentation'  errors),  whereas,  Gu,  reflects  only  the  geometry  of 
the  user  systems  and  is  the  nx4  matrix  of  partials  where  n  is 
the  number  of  visible  satellites.  [Milliken  and  Zoller  I960}. 

The  process  for  determining  this  purely  geometric  effect  on  a 
GPS  user  position  in  each  axis  is  described  below  step  by  step. 

(1)  Change  the  station's  geodetic  coordinates  0,  A  ,h,  to 
earth-fixed,  earth-centered  cartesian  coordinates,  X,  Y,  Z. 

X  »  (N  +  h)  cos  ^  cos  A 

Y  -  (H  +  h)  cos  <P  cos  A 
Z  »  ((l-e2)N  +  h)  sin.  ,  where: 

N  is  spheroidal  normal  length  terminated  by  earth's  minor  axis. 

N  *  a/(l-e2sin2  ♦  )l/2, 

stnd  a  is  the  semi-major  axis  and  e  the  first  eccentricity  of  the 
reference  ellipsoid. 

(2)  Calculate  the  direction  normal  to  ellipsoid 
N  is  the  unit  vector  in  the  direction  normal  to  ellipsoid 


» 

• 

cos4  cosA 

"i 

N  - 

cos 4  sinA 

9 

n2 

sin* 

^  ■ 

t?3 

►  .4 
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(3)  Calculate  the  topocentric  range  vector 

given *  Xata,  Yata,  Zata  are  the  earth-centered,  earth- fixed 
cartesian  coordinates  of  the  station. 

Xsv  ,  Yav,  Zav  are  the  earth-centered,  earth- fixed 
cartesian  coordinates  of  the  satellite. 

then:  Tt  is  the  range  vector 

^sv  ~  ^sta 
Y#v  ~  *sta 
zsv  “  zsta 

(4)  Calculate  the  range 
given:  Tt  is  the  range  vector 

then:  |r|  is  the  range  (distance  from  station  to  satellite) 

I  a|  -  (if  -KR*  +  R^)1/2 

I  al  -  ((I„  -  -  w2  *  (2„  -  z,^bui 

(5)  Find  the  direction  cosine  (unit  vector)  of  the  range 
vector 

given:  R  is  the  range  vector 
]r|  is  the  range 

then:  TT  is  the  direction  cosine  (unit  vector)  of  the  range 
vector. 


(6)  Find  the  azimuth  of  the  satellite 

given:  K  the  rotation  matrix  to  convert  XYZ 

coordinates.  -cosA  slat  -sinAsint 

*  •  -siaA  cosA 

cosAcost  slnAcos t 


U  is  the  direction,  cosine  of  the  normal  to  the  reference 


ellipsoid. 

To 
*3 

then:  az  is  the  azimuth  of  the  satellite 
an  *  arctan  (T2/T1) 

(7)  Calculate  the  elevation  angle  of  the  satellite 
given:  U  and  ft 

find:  *-T  m  +>  ^2^2  + 

and  V-f  •  (fll-ltj  cosd  ■  coed 

then:  EA.  is  the  elevation  angle  of  the  satellite 
EX  *  90°  —  9r  9  *  90°  —  E2U 


find: 

a?)  cu)  - 


Figure  2-3— L  Geometry  of  Elevation  Angle 


coa ( 90°  -  EA)  -  U]Ni  +  U2»2  +  U3N3 
sin(EA)  -  UjNi  +  U2N2  +  U3N3 
Therefore,  EA  «  ARCSIN  (U^Ni  +  U2N2  +  U3N3) 

(9)  Solve  for  Gu  matrix 

given:  Un  is  the  direction  cosine  of  the  range  vector  of 
satellite  n. 

b  *  cat  is  the  range  equivalent  of  user  clock  offset  (At) 
and  n  is  the  number  of  visible  satellites. 

.  C  is  equal  to  the  speed  of  light, 
then:  Gu  is  the  n  X  4  matrix  of  partials 


aCu4 
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where  Rn  *  C(xn-xatn)2 

3*u  , 


(2.3  -3) 

(Yn-ystn)2  +•  (2n2-*stn)23x/2  +■  b 


and  therefore 


1. 


The 


is  the  observation  error 


equal  to  1  meter. 


where: 


(9)  Compute  Covariance  matrix,  COV 


°o2  (sJCu)'1 


Cov 


2  T  -1 
°o(G  G  ) 
u  u 
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2 

XX 

cxy 

c 

xz 
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y* 
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zx 
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cby 
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2 

where  c.  .  is  equal  eo  aa<*  cli"aii’ 

.  .  ,  2 

3  ii  *»**  variance  of  the  range  obaervation. 

(10)  Calculate  POOP  from 
(G^Gu)”1  the  4X4  matrix 
PDOP  *  (  xx2  +■  yy2  +■  zz)l/2 

(11)  Find  GDOP 

given  (G^GU)  1  is  the  4X4  matrix 

l/  2 

then:  GDOP  -  (<r2  «*►  <j 2  «♦«  c2,  a!.)  meters 

xx  yy  zz  bb 

2.4  Station  Tracking  Geometry 

The  first  important  graph  (figure  2.4-2)  of  Elevation  angle 
versus  time  of  each  GPS  vehicle  used  in  this  study  enables  the 
user  to  know  the  best  time  to  track  these  satellites  from'  a 
ground  station  anywhere  in  the  world.  This  also  verifies  the 
GDOP  values  versus  time  by  empirical  observation. 

These  Azimuth  versus  Elevation  Angle  plots  verify  the  GDOP 
and  PDOP  calculations  by  showing  the  two  viewing  periods  during 
a  24  hour  period  in  a  "bulls  eye"  plot  pattern  (see  Figure 
2.4-1).  This,  therefore,  enables  the  user  to  select  the  optimal 
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FIGURE  2.4-1 
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3.  VARIANCE-COVARIANCE  STUDY 
3. 1  Introduction 

In  this  chapter,  the  development  of  the  least  square* 
estimation  theory  and  data  acquisition  for  batch  processing  is 
described.  In  Section  3.2.,  the  definitions  and  justifications 
for  both  the  mathematical  models  and  parameter  set  are  pre— 
sented.  The  technique  of  estimation  for  least  squares  adjustment 
and  the  procedure  for  linearizing  the  mathematical  model  are 
described  in  Section  3.3.  Section  3.4  illustrates  the  procedure 
that  utilizes  the  consider  covariance,  that  is,  the  method  for 
analyzing  the  effects  of  unestimated  parameters.  Finally, 
Section  3.5  describes  apriorl  covariance  statistics  and  the 
partial  derivatives  used  in  the  process  to  establish  the 
covariance  analysis. 

The  adjustment  philosophy,  [Bierman,  1977]  was  developed 
at  the  Jet  Propulsion  Laboratory  (JPL)  and  was  incorporated  in 
their  tracking  and  orbit  determination  section.  The  choice  of 
the  parameters  set  in  our  study  was  influenced  greatly  by  past 
studies.  [e.g.  Fell,  1980].  Therefore,  the  emphasis  in  our  study 
is  on  the  selection  of  unestimated  parameters, i.e.  consider  pa¬ 
rameters  and  their  influence  on  the  estimate  parameters.  The 
computation  of  the  consider  covariance  enables  the  user  to 
estimate  the  errors  caused  by  incorrect  apriori  statistics  or 
mismodeling.  The  intention  is  to  optimally  analyze  the  effects 
of  the  unestimated  or  consider  parameters  so  that  the  propaga¬ 
tion  of  error  covariances  in  the  estimate  parameters  are  readily 
available  and  can  be  easily  computed.  An  important  feature  of 


estimation  algorithms  is  that  error  covariances  for  the  esti¬ 
mates  are  readily  available  or  can  be  readily  computed.  This 
covariance  knowledge  is  used  to  measure  confidence  in  the 
estimate  parameters  and  discrepancies  between  them.  The 
consider  covariance  arises  in  the  evaluation  of  model  that  omits 
certain  parameters,  and  we  speak  of  considering  the  effects  of 
these  unestimated  parameters,  thus  called  consider  parameters. 
CBierman,  1977] 

3. 2  Development  of  The  Mathematical  Model 

This  section  describes  the  mathematical  model  for  the  least 
square  adjustment  using  estimation  theory.  The  mathematical 
model  is  suited  for  GPS  pseudo-range  measurements  and  the 
optimal  utilization  of  estimate-consider  covariance  analyses. 
The  mathematical  model  is  developed  from  previous  studies  CFell, 
1980],  and  is  implemented  so  that  it  may  be  easily  expanded  to 
meet  requirements  for  future  use.  The  choice  of  a  parameter  set 
is  justified  by  conclusions  from  previous  studies  about  GPS 
pseudo-range  measurements . 

3.2.1  Mathematical  Model 

A  GPS  observer  vat  ion  for  geodetic  point  positioning 
requires  pseudo  range  measurements  from  four  GPS  Space  Vehicles 
with  time  being-  the  fourth  solution  variable. 


Let  range  Ri  be. 


aL  *  ( <usi T  ttstaP  2  (,r3iT  Tsta^ 2  +  "sta^  2) 1/2  (3.2-1) 

the  topocentric  range  from  any  ground  station  to  any  satellite 
position  is  shown  in  Figure  3.2.1,  where  the  earth  fixed 
coordinate  system  (u,v,w)  is  oriented  towards  Greenwich  mean 
astronomical  meridian  (u  axis)  and  the  Conventional  Internation¬ 
al  Origin  (w  axis)  with  v-axis  forming  a  right-handed  coordinate 
system  with  u  and  w,  this  coordinate  system  being  defined  by 
Bureau  International  de  1‘Heure  (BIH)  [Fell,  1980,  page  30]. 


The  accuracy  of  satellite  ephemerides  and  tropospheric 
refraction  modeling  and  the  stability  of  satellite  and  receiver 
clocks  will  have  important  consequences  in  the  application  of 
range  observations  to  geodetic  positioning.  [Fell,  1980,  Pg. 
22].  From  this  type  of  statement  the  mathematical  model  adopted 
for  this  covariance  analysis  is  formulated  to  result  in  the  fol¬ 
lowing  equation.  Opr,  the  theoretical  pseudo  range  observation, 

0pr  *  o  +•  a*trop  +>  c  (b  +  dr*t  +  ag*t2)  (3.2  -2) 

where  p  is  the  computed  topocentric  range 

a  is  scaling  factor  of  tropospheric  correction 
trop  is  the  actual  tropospheric  correction 
c  *  value  for  speed  of  light 

*  t-tQ  where  t0  is  the  time  of  first;  observation  and 
in  this  study  tQ  »  0. 
b  »  combined  satellite-station  clock  bias 
dr  *  combined  satellite-station  clock  drift,  first 
derivative  of  clock  error. 


ag-  *  combined  satellite-station  clock  aging,  second 
derivative  of  clock  error. 

3.3  Least  Squares  Estimation  Theory. 


The 

utilization  of  the 

estimation 

theory 

for 

least 

squares 

adjustment  begins  with 

a  review 

of  the 

"classical" 

linear 

least  squares  problem. 

[Uotila, 

1967]. 

Next, 

the 

linearization  of  a  non-linear  mathematical  model  is  described. 


Also,  a  detailed  description  o£  the  partial  derivatives  for  the 
Jacobian  matrix  or  the  "coefficient  matrix  of  observation 
partials"  [Bierman,  19773  is  included. 

Suppose  we  have  a  linear  system, 

z  *  Ax  +  V  (3.3-1) 

where  z  is  an  m  vector  of  observation,  X  is  the  n  vector  of 
variables  that  are  to  be  estimated,  A  (m,n)  is  the  'coefficient 
matrix  of  observation  partials  and  V  is  an  m  vector  of  observa¬ 
tion  errors. 

The  least  squares  solution,  X^g,  to  the  linear  system  is 
chosen  to  minimize  the  mean  square  observation  error:  ■* 

J  »  ,£  v2  -  vTv  (3.3  -2) 

3-1  3  . 

J(x)  »  ( r  —  Ay)T(z  -  Ax)  (3.3  -3) 

i .  e . , 

When  more  than  one  X  minimizes  the  mean  square  error  J,  the 
X^s  is.  chosen  to  minimize  X  of  smallest  Euclidean  length  denoted 
as, 

H*  (3-3-4) 

and  this  uniquely  defines  X^g. 

Therefore,  J(x)  is  nonnegative  and  quadratic  in  the  compo¬ 
nents  of  X  so  that  a  necessary  condition  for  a  minimum  is  that 
the  first  variation, 5 J(x) ,  vanishes. 
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The  first  variation  is  defined  as, 

J(x  +  <Sx)  -  J(x)  -  5 J (x)  +  o(^5x)  (3.3  -5) 

where 

o(5x)  is  such  that 

lim  (o  ( 5x)  /  |j  5x|| )  -0 

$x-0  (3.3-6) 

5  «J(x)  is  obtained  by  applying  rules  of  differentiation  and 
matrix  manipulation,  thus, 

ijc*)  -  S[<z  -  Ai)I](r  -  Ax)  +•  Cr  -  Ax)TJf*  -  A,] 

"  (-Aix)T(r  -  A_)  +•  (z  -  A_)T(-A4x) 

*  X  ' 

•  5xT(ATAx  -  ATz)  +  r5xT(ATAx  -  ATz)]T  (3.3-7) 
For  the  first  variation  &J(x)  to  vanish  for  all  varia¬ 
tions  fix,  it  is  necessary  that  x  satisfy  the  normal  equations, 

ATAX  •  ATz  (3.3  -0) 

It  is  important  to  note  when  m  n  and  A  has  rank  n,  then 
AtA  is  non-singular  and.  the  solution,  Xls ,  is  unique,  i.e., 

Xls  *  (ATA)— LaTz  (3.3  -9) 

Mow,  suppose  we  are  given,  a  non-linear  system  in  the  form 

r  -  *  (3.3-10) 

with  nominal  values  assigned  to  the  estimate  and  consider  Xe,  X^, 
parameters,  where  Y  is  an  m  vector  ot  observations,  and  X  is  the 
n  vector  of  variables  tha-t  are  to  be  estimated  or  considered. 

The  linearization  of  the  mathematical  model  is  obtained  in  a 
Taylor  series  expansion  such  that: 

*(**•*<.)  ~  +  3£_  +■  H . 0 . T.(3. 3—11 ) 

(H.O.T.  means  higher  order  terms)1-  J 


Writing, 


1 


we  have, 

y  -  f # (x£ , x£) 


+  H.O.T.  +  v 
3x 


(3.3  -12) 


Considering  that  XP  is  close  to  X  the  higher  order  terms  may 
be  neglected,  thus, 


Xf  •  Y  -  f  • 


3  £ 
33T 


At  +  ~v 


(3.3-13) 


Mow,  using  the  notation  of  y  for  Ay  « 


and  x  for  ,  i.e.. 


:.Tf 

- 


We  may  write: 

t  -  kUJ 

a 

*e 

L  ••  3 

+  v“ 


(3.3  14) 


or 

r  •  (3.3-15) 

the  outcome  of  a  non— linear  model  for  estimation  theory  of  least 
squares  adjustment. 

3.4  Development  of  Consider  Covariance. 


In  this  section  is  presented  the  derivation  of  the  consider 
covariance  that  arises  when  we  speak  of  considering  the  effects 
of  the  unestimated  parameters, consider  parameters,  X^.  The 


estimate  parameter a ,  and  the*  truth  nodel  being  defined  as 

T -  A^re  +  A^c  +  "v  (3.3-16) 

with  the  assumption  that  Ae,  Ac,  Xe  are  deterministic.  There¬ 
fore, 

E(xe)  *  Xe  *  a  constant. 

Also,  assume  that  Xc  and  V  are  random  variables  with  known 
covariance  and  the  equations  have  been  nomalized  to  show  the 
following  characteristics  so  that  * 

E(V)  *  0  and  E(WT)  -  Zv  (3.3-17) 

where  Zv  is  the  m  dimensional  identity  covariance  on  measure¬ 
ments. 

Let  EGrJ  » jr  •  ,  where  is-  the  average  value  of  the 

variable  over  all  possible  values. 

Then,  ElCic -irc)(3^ -irc)r]  -Tt  (3.3-18) 

where  £c  is.  the  coneidec  covariance. 

By-  simple  manipulation-  of.  3~l-ia  we  gets 

xctye^j  •  rc  +.  (3.3-19) 

We  assume  that  K(Tzf)  —  0  •  that  is,  there  is  no  correlation 

c. 

between  the  errors  on  observations  and  the  consider  parameters. 
Let:  the  filter  model  be, 

^  (3*3-20) 

and  frost  the  derivation  of  the  least  squares,  solution  we  have, 

'  5*.  -  [a$vlv'la£;V  (3*3-2l) 

-  (Vf.  *  *> 

*  K4*1*  A«^  <4**i.  +■  *  v) 

"  t4^*1’4J  L*«S»-4«**  ■*'  ■*■••• 


(3.3-22) 


If  we  now  define  the  following  matrices, 


'  [A^f1 

s,  -  P^r;1 

then  (3. 3. -17)  reduces  to: 

xe  *  +  ScXc  +>  S*v 

**  "  **  *  SjJ*,.  -K  svv 

Examining'  the  statistical  characteristics  of  random  variable  X; 
we  have 

Et^  **  *e]  -  E(S«>  -t%-  S^d,.)  +  g^zc?) 
from  the  fact  that  E(c)  »  c  where  c  is  a  constant* 

If  Xc  is  biased,  i.e.,  E  (Xc) *0,  then  Xe  is  a  biased  estimate 


(3.3-23) 

(3.3-24) 

(3.3-25) 

(3.3-26) 


therefore,  E(x%)  -  £ 


Since  the  estimate  covariance  may  be  defined  as. 


2*e  *  E[(*e  -  *e)(Xe  -  f*)7]  using  equation  (3.3-26) 

•  E[ (ScSc.  +  M><Mc  4-  §V?)T] 

-  E[ScxcsJ§c  +  Scic^S?  4-  S>**S£  4-  S>*T§*] 

-  §aE(2y£)§*  4-  ScE^T)§J  4-  S*E(fc£)§*  +  SVE(^T)SJ 

-  §£  <rc  +  nyf)^  +  uain9  < 3-  3~19  > 


+  scuc^sc  4.  v:  V 

scrcs£  4-  scucl£§£  4-  Pep;i?a 


Tr_1v  .-1*  S 


e  a 


and  (3.3-15) 


using  equation  (3.3-25) 


therefore#  the  estimate  covariance  may  be  written  in  the  form# 
22e  *  4-  Scu^uJ§T  4-  (3.3-24) 

and  Sc  defined1  in  equation  (3.3—24) 

If  Xc  La  unbiased,  i.e-,.  u*c  «  & 

Therefore#  the  total  covariance  in  finally  derived  to  include 
the  consider  covariance# 

IS.  -  ?.4-S^r tSl  (3.3-2S) 


3.5  Partial  Derivatives: 

Sow-  that  the  covariances  have  been  derived#  we  need  to 
consider  the  partial  derivatives  for  the  Jacobian,  matrix.  The 
interactive  program  is  capable  of  handling-  multiple-  variables 
from  the  mathematical  model  for  pseudo  range  of  measurements. 
This  section  describes  the  technique  for  solving  the  adjustment 
algorithm  involved  in  covariance  analysis.  In  the  following. 


sections,  the  discussion  will  contain  information  about  certain 
weight  matrices,  apriori  statistics,  and  partial  derivatives 
involved  in  the  method  for  the  estimation  theory. 

Assume  Opr,  theoritical  pseudo  range  measurement,  see 


equation  (3. 2.-1)  to  be  given  by  the  mathematical  model. 


Opr  *■  o  +■  a*trop  +  c  (b  +-  dr*t  +  ag*t^) 
where  the  variables  are  defined  in  Section  3.2. 


(3.5-1) 


The  list  of  parameters  for  this  particular  study  is  as  follows: 
NSTA  latitude  of  station  (meters) 

ESTA  longitude  of  station  (m) 

USTA  height  of  station  (m) 

TROP  tropospheric  correction  (m) 

LSAT  long-track  of  satellite  (tangential)  (m) 

CSAT  cross-track  of  satellite  (m) 

RSAT  radial-track  of  satellite  (m) 

BIAS  combined  user  and  satellite  clock  bias  (Sec) 

DRIFT  combined  user  and  satellite  clock  drift  (Sec/ sec) 

AGEING  combined  user  and  satellite  clock  aging  (Sec/sec2) 

This  study  uses  the  NSTA,  ESTA,  USTA  as  estimate  parameters 
and  shows  the  correlation  between  any  of  the  other  parameters 
and  their  effect  on  the  NSTA,  ESTA,  and  USTA  for  geodetic 
positioning  techniques  specifically  pseudo  range. The  impact  of 
consider  parameters  onto  the  estimable  parameters  and  the 
numerical  degradation  (failure)  in  estimate  correlated  parame- 
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3.5.1  Partial  Derivatives  Development. 

He  first  consider  the  measurement  model  involving 
pseudo-range/  PR/ 

0pr  ■  p.  +  a*trop  +•  c  (b  +■  dr*t  +*  ag*t^) 


PR 


The  partials  are  computed  in  a  subroutine  PDERIV  and  stored  in 
an  array  for  each  iteration  in  the  following  fashion.  (Chain 


rule) 


3°pi 


3NSTA 

30pr 


3ESTA 

30pr 


30STA 


•  _ 

• 

3  XSTA 

3XSTA 

3NSTA,E5TA,USTA 

** 

3  3 

m 

(3.5-2) 


where  Xsta  the  position  vector  of  the  station,  these  partials 
result  in  being  equal  to 


(3.5-3) 


where  line  of  sight,  L.O.  S. ,  is  equal  to  a  unit  vector,  or 
sometimes  referred  to  as  the  direction  cosine,  and  is  defined  as 
a  normalized  vector, 

L.O.S.  »  Syr  -  Sacs  (3.5-4) 

length  of  (Rsv-Rsca) 

where  RSV  and  RSTA  are  the  position  vectors  of  satellite  and 
station,  respectively.  ROTSTA  is  a  rotation  matrix  transforming 
North,  East,  and  Up  coordinates  (local  station  coordinate 
systems  see  Figure  3.5-lb)  to  earth-fixed  cartesian. 
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The  next  partial  of  the  mathematical  model ,  ,  that  is,  the 
partial  of  Opr  w.r.t  alpha, a,  the  scale  factor  of  the  tropos¬ 
pheric  correction  is,  30nr 

p  -  *  TROP 

(3.5-5) 


3a 


where  TROP  is  the  computed  tropospheric  correction  and  is  found 
using  the  CHAO  model. (1983)  The  CHAO  model  (Reference,  1983)  was 
tested  and  evaluated  by  the  GST  department  at  DMA. 

The  partials  for  long-track,  cross-track  and  radial-track 
consider  parameters  aura  computed  in  a  similar  method  as  the 
ground  station  ones.  Thus,  we  have. 


30pr 

3LSAT 

30pr 

3CSAT 

30pr 

'  M  ‘  3 

3XSAT 

3LSAT,CSAT,RSAI 

w 

3  (3.5-6) 

3BSAT 

4 

where  SATROT  * 


(3.5-7) 


and  1,  c,  r  are  the  components  for  a  local  satellite  coordinate 
system,  (see  figure  3. 5 -la) 


7XGUBE  3  -5— 1*» 

Local  eafllita  coordinate,  sys 


The  process  for  the  computation  is  as  follows,  where  R  * 
that  is ,*R, the  unit  vector  of  the  satellite,  where  R*  (x,y,z) 
of  the  satellite. 

/>  R 

Next,  set  up  a  vector, R  “-vrwhere  r  «  (x,y,z)of  the  satellite. 

|R| 

that  is,  the  unit  velocity  vector  of  the  satellite. 

A  A 

The  process  continues  by  producing  the  cross  product  of  R  X  R. 


A  ♦ 

which  is  c  «  R  x  R 


Finally, 

1  *  c  X  r  (cross  product) 
Thus  we  have. 


3°pr 

T 

,  , 

3LSAT 

1  1 

p  T 

-  L.O.S .  *• 

1  ^  c  r 

3C5AT 

lL  J3 

3<3pr 

1  i 

1  a 

3BSAX 

3 

(  1 

The  parti als  for  the 
found.  Therefore, 


?°2S _ _ 

3BIAS 


bias,  drift,  and  aging  terms  are  easily 


3°pr  .  c  30pr 

-  az»iy  *“  '3agbg 


e* 


where  t  represents  the  time  of  the  pseudo  range  measurement. 
3.6.2  Apriorl  covariance  matrices  and  apriori  uncertainties  for 
measurement  types  and  parameters. 

The  weight  matrix,  called  COW,  is  set  up  for  the  three 
measurement  types  pseudo  range,  range  differencing,  and  second 
range  differencing  is  as  follows: 


r . 


T 


(3.5-9) 


where  la  =*  NMEAS  *  1  the  pseudo  range  measurement  type  for  this 
study.  The  interactive  program  requires  the  number  of  measure¬ 
ment  types  for  the  adjustment  algorithm  and  may  account  for  all 
three  measurement  types  in  the  future.  COVE0  is  the  apriori 
covariance  matrix  for  estimate  parameters  and  is  formed  in  this 
fashion. 


COVE0 


r 


(3.5-10) 


where  j  *  1,  .  .  .  ,n  the  number  of  estimate  parameters. 

COVC  is  the  apriori  covariance  matrix  of  consider  parameters, 
and  is  formed  this  way. 


(3.5-11) 


where  i  =  1,  .  .  . 
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n  the  number  of  consider  parameters. 


The  a  priori  uncertainties  for  the  measurement  types  and 
parameters  were  found  by  investigating  recent  studies  (personal 
communication)  performed  by  Naval  Surface  Weapons  Center,  NSWC, 
Dahlgren,  Virginia.  There  is  a  table  of  these  apriori  values 
found  in  the  appendix. 


4.  EXPERIMENTAL  RESULTS 


4.1.  Test  Plan  Design 

The  experiments  involved  in  this  study  show  the  effects 
of  each  consider  parameter  individually  onto  the  estimate 
parameters  by  means  of  a  consider  covariance.  The  experimental 
test  plan  for  this  study  involved  the  performance  of  fifteen 
experiments  with  results  displayed  on  graphs.  The  least  squares 
adjustment  algorithm  allows  the  user  to  chose  estimate  and 
consider  parameters.  The  first  seven  experiments  in  this  study 
use  the  tracking  station  Austin,  Texas  and  the  estimate  parame¬ 
ter  are  NSTA,  ESTA,  and  USTA  which  in  the  E-N-U  coordinate 
system  is  C Conley,  19842  is  adopted  for  this  study.  The 
experiments  use  the  five  GPS  vehicles  currently  operational  with 
the  pseudo  random  noise  (PRN)  numbers  4, 5,6,8  &  9.  The  epoch  cf 
the  orbits  is  the  same  for  all  orbits  Day  17,  0  seconds,  1983. 

The  first  experiment  involves  the  solution  set  of  the 
latitude4  (N),  longitude  (E)  and  vertical  (U)  components  of  the 
local  coordinate  system  for  the  station  Austin,  Texas.  In  Table 
1,  experiment  number  1  the  results  show  the  sigmas  of  latitude, 
longitude,  and  vertical  being-  .1,  .2  and  .1  meters  respectfully* 

The  second  experiment  calls  for  the  solution  set  of 
estimate  parameters  NSTA,.  ESTA  &  USTA  and  shows  the  effects  of 
the  consider  parameter,  TROP,  onto  the  estimate  parameters, 
tropospheric  correction  by  means  of  the  consider  covariance 
(formula  3.3-1).  This  experiment  shows  what  happens  to  each 


estimate  parameter  when  the  tropospheric  correction  is  consid¬ 
ered,  thus  the  name  consider  parameter.  Each  experiment  solves 
for  NSTA,  ESTA  and  USTA  and  uses  the  remaining  parameters  in  the 
model  such  as  ephemerides  parameters  (LSAT,  RSAT,  CSAT)  and 
combined  satellite-station  clock  parameters  (BIAS,  DRIF,  AGNG ) 
plus  the  tropospheric  correction  individually  as  consider 
parameters.  The  model  of  pseudo-range  observations  from 
tracking  satellite  to  station  is  the  one  developed  in  this  study 
as  formula  3.2-3.  The  observation  error  used  for  pseudo  range 
is  assumed  to  be  1  meter.  The  assumed  standard  deviations  used 
in  the  experimental  tests  from  the  list  of  parameters  on  page  30 
are  as  follows: 


NSTA 

XX 

50  meters 

ESTA 

« 

50  meters 

USTA 

m 

50  meters 

LSAT 

as 

8  meters 

RSAT 

« 

2 -meters 

CSAT 

m. 

4-  meters 

BIAS 

m 

10“5  gee. 

DRIF 

» 

10*10-  sac /sec 

AGNG 

10-18  sec/sec 

The  test  results  from  the  experiments  indicate  the 
reliability  of  GPS  positioning  with  respect,  to  changes  in  the 
geometry  with  time  and  the  effects  of  consider  parameters.  The 
changes  in  geometry  were  indirectly  measured  by  the  changes  in 
the  standard  deviation  of  the  estimate  parameters  and  changes  in 
the  correlation  coefficients  among  the  estimate  parameters. 


To  summarize,  the  following  factors  stipulate  each 


experiment : 

♦Constants 

-five  GPS  satellites 

-tracking  station  Austin,  Texas 

-Pseudo-range  observation  error,  1  meter. 

-Elevation  angle  cutoff  15  degrees. 

♦Estimate  parameters  in  the  experiments  ESTA,  NSTA,  USTA, 
♦Consider  or  unestimated  parameters  are  satellite  long,  cross 
and  radial  track,  and  tropospheric  correction,  and  clock 
parameters  bias,  drift,  and  aging. 

4.2  Test  Plan  Sequence  for  Each  Experiment. 

Preprocessor'  activity  for  each  experiment  begins  with  the 
position,  velocity,  and  ground  trace  data  for  each  satellite. 
The  azimuth  and  elevation  angles  versus  tima  of  the  range 
observations  from  each  satellite:  are  computed  and  stored  in 
files  using  the  UNIVAC  1100  series  computer.  This  data  is  used, 
in  the  station  geometry  study  and  found  to  be  useful  in  the 
evaluation  and  analysis  of  each  experiment.  The  files  are  then 
transformed  into  plotable  files  for  the  Hewlett  Packard  HP  9845 
desktop  computer.  The  output  display  is  found  in  Section  2.2 
page  .The  POOP  and  GDOP  plots  as  a  function  of  time  show  that 
the  number  of  satellites  in  view  and  the  tracking  time  for  the 
receiver  is  sufficient  for  a  solution.  The  GDOP  study  shows  the 
user  that  the  optimal  time  for  satellite  tracking  is  from  16d 
12h  to  16d  16&. 
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4.3  Processor  Activity  £or  Each  Experiment. 

The  parameters  NSTA,  ESTA  and  USTA  remain  estimate 
parameters  for  each  experiment  and  the  experiments  are  comprised 
of  the  consider  parameter  set  using  one  at  a  time  making  seven 
experiments.  For  instance,  the  second  experiment  shows  the 
effect  of  tropospheric  correction  (TROP)  onto  the  estimate 
parameters  NSTA,  ESTA,  and  USTA. 

The  second  experiment  results  found  in  Table  1  shows  a 
0.1  to  0.3  meter  difference  in  the  standard  deviation  for  the 
latitude  component  NSTA.  The  standard  deviation  of  ESTA  shows  a 
.3  to  .8  meter  difference  caused  by  the  consider  parameter, 
tropospheric  correction.  The  USTA  standard  deviation  took  the  a 
change  being  from  .4  meter  to  .5  meters. 

The  third  experimenton  Table  1  is  the  NSTA,  ESTA,  and 
USTA  as  estimate  parameters  and  long^  track  of  the  satellite 
(LSAT)  as  a  consider  parameter.  In  the  final  solution  there  is 
1.2  meter  difference  shown  on  the  NSTA  parameter.  The  same 
result  is  found,  in  the  ESTA  parameter.  The  major  difference 
found  in  the  USTA  parameter  is  .8  meter  caused  by  the  consider 
parameter ,  LSAT . 

The  fourth  experiment  on  Table  1  is  the  NSTA,  ESTA,  USTA 
as  estimate  and  cross-track  of  satellite  (CSAT)  being  the 
consider  parameter.  The  consider  parameter,  CSAT,  had  an  affect 
on  the  final  solution  for  NSTA  from  1.6  to  2.2  meters.  The 
consider  parameter  showed  approximately'  a  .  5  meter  noise  level 
in  the  estimate  parameter  ESTA.  The  consider  parameter  had  a 
negligible  effect  on  the  USTA  estimate  parameter. 


The  fifth  experiment  on  Table  2  involving  the  consider 
parameter,  radial-track  of  satellite,  (RSAT)  affected  the 
estimate  parameter  by  approximately  .2  to  .8  meter.  The  ESTA 
parameter  took  a  1.0  to  2.2  meter  jump  from  the  RSAT  consider 
parameter.  The  RSAT  has  a  0.4  to  1.8  affect  on  the  USTA 
parameter. 

The  next  experiments  involving  the  BIAS,  DRIFT  and  AGING 
consider  parameters  revealed  the  most  significant  influence  on 
the  state  parameters. 

The  fifth  experiment  involves  the  station-satellite  BIAS 
clock  consider  parameter.  The  most  significant  errors  are  found 
in  the  state  parameters  from  this  consider  parameter.  All  three 
estimate  parameters  show  a  degradation  but  the  USTA  parameter 
took  the  most  "beef"  from  the  consider  parameter. 

The  sixth,  experiment,  involves  theconsider  parameter, 
clock  drift,  DRIF,  in  which  a  10  meter  fluctuation  in  the  state 
parameter  N5TA.  The  ESTA  parameter  shows  a  10m  fluctuation  due 
to  the  consider  parameter.  The  USTA  has  the  most  significant 
fluctuation  with  a  600  meter  difference. 

In  the  seventh  experiment  the  station-satellite  clock 
aging  shows-  negligible  fluctuation  from  the  state  parameters. 
The  only  parameter  that  shows  a  .1  meter  fluctuation  due  to  this 
consider  parameter  is  the  USTA  parameter  and  a  station  tracking 
geometry  useful  for  simulation  studies.  The  point  is  made  that 
there  is  a  level. of  confidence  added  to  the  measures  of  GDOP.  At 
this  level  the  study  is  very  useful  for  setting  up  tracking 


station*  testing  receiver  capabilitiss.  Tha  chaptsr  is 
written  in  a  fashion  so  that  an  inexperienced  sciantist  may 
undar stand  tha  procass  for  computing  azimuth  and  alevation 
angles  of  a  satallita  from  a  station  on  tha  earth's  surface. 
Also#  there  is  a  short  description  on  tha  GEOSTAR  receiver 
currently  being  tasted  at  tha  University  of.  Texas  by  Applied 
Research  laboratory  (ARI»)  at  Austin#  Texas. 

4.4  Correlation  Coefficients 

Experiments  nine  through,  fifteen  found  on  Tables  3  and  4 
shoe  the  interrelationships  among  the  parameters  involved  in  the 
math  model  for  pseudo  range.  A.  value  of  unity  indicates  that 
the  parameters  are  highly  dependent  and  therefore  produces  an 
unreliable  system  for  improvements  in  station  coordinates  for 
point-positioning  using  the  pseudo  range  technique. 

This  study  found  a  high  correlation  among  certain 
parameters  such  as  the  tropospheric  correction  has  &  .8  correla¬ 
tion  coefficient  with  the  vertical  component  (Table  3).  Another 
significant  finding  is  a  positive  correlation  coefficient  with 
cross-trach  of  satellite  and  the  longitude  and  vertical  compo¬ 
nents  of  the  station..  Also  in  Table  3  the  radial-trade  of 
satellite  and  the  vertical  component  of  the  station  have  a  .8' 
correlation  coefficient.  Finally#  in  Table  4  there  is  a  high 
correlation  coefficient  between  bias#  drift  dock  parameters  and 
the  vertical  component  of  the  station. 


5.  CONCLUSIONS  AND  FUTURE  RESEARCH 
5.1  Summary  and  Conclusions 

A  GPS  Pseudo-range  covariance  analysis  of  the  interactive 
program  is  presented  in  Appendix  A  for  station-tracking  geome¬ 
try,  experiment  planning,  and  future  design  analysis.  Tables  and 
figures  for  ease  of  operation  and  adaptation  are  presented.  The 
program  itself  is  well  documented  for  future  users  (example; 
expanding  parameter  set).  By  an  explanation  of  the  theory  on 
which  the  program  is  based  and  of  the  mathematical  models  Which 
it  incorporates,  an  overview  of  GPS  pseudo-range  process  is 


given. 


An  introductory  chapter  touches  upon  the  historical  and 


future  aspects  of  GPS  as  well  as  its  application,  to  point-posi¬ 
tioning. 

Chapter  2  opens  with  the  process  for  determining  the 
preprocessor  stages  that  is  the1  generation  of  the  GPS  satellite 
ephemeris,  the  subsatellite  points  projected  on  a  world  map, 
azimuth  and  elevation-angle  plots  of  the  satellite  from  ground 
station  versus  time,  and  the  GDOP  study.  A  description  of  the 
method  for  computing  azimuth  and  elevation  angle  of  a  satellite 
from  a  ground  station  is  presented  completely  in  a  tutorial 
fashion.  The  development  of  the  geodetic  receiver  (GEOSTAR)  is 
presented  along  with  a  broad  description  of  its  capabilities. 

In  Chapter  3,  the  development  of  the  least  squares 
estimation  and  the  method  used  for  the  covariance  analysis  is 
described.  The  process  and  mathematical  models  are  described 
for  the  development  of  consider  covariance  and  the  partial 


derivatives.  Clearly  more  research  is  needed  to  improve  the 


mathematical  models  and  the  estimation  process.  In  the  next 
section,  an  attempt  is  made  to  qualify  the  future  areas  of 
research. 

5.2  Suggestions  for  future  research. 

This  research  will  no  doubt  continue  since  the  need  for 
accurate  and  timely  data  analysis  is  accelerating.  several 
suggestions  that  may  help  to  enhance  the  future  needs  are:  l)to 
test  different  elevation  angle  cutoffs  for  the  best  observations 
and  2) to  experiment  with  the  tropospheric  corrections  models. 
Additionally,  there  could  be  an  attempt  to  determine  the 
relationship  between  time  bias  and  the  solution  for  height  of  a 
station  along  with  different  elevation  angle  cutoffs-  One  could 
study  the  effect  of  3rd  order  ionospheric  errors  and  how  they 
interact  with  the  estimate  parameters.  Also,  a  question  arises 
as  to  the  advantages  pseudo- range  measurements  have  over  Doppler 
or  range  differencing  methods.  As  more  satellites  are  launched, 
there  should  be  an  attempt  to  have  a  selective  method  to  obtain 
the  "beat"*  satellite  observations  from  passive  tracking  sta¬ 
tions.  Further  consideration  should  be  given.  to  the  cost 
effectiveness  of  these  proposals  for  GPS  point  positioning. 


.V.v. 


56 


A-.' 
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ESTIMATE  CONSIDER 

parameter  parameter 


0.1«j<0.3  0.3<a<0.8  0.4<a<0.5 


1.  2«7<L.  8 

1.  6<cr<2.  9 

1.0«X<1*8 

1.6<<7<2.2 

0.  4<0<2.  0 

1. 4«y<Z.  0 

S  Latitude  Component 
S  Longitude  Component 
U  Vertical .  Component 
T  Tropospheric  correction 
L  Long-track  of  satellite 
C  Cross-track  of  satellite 
R  Radial- track  o'f  satellite 

(  <7  units  in  meters) 


TART.R  2 


Variations  in  o  when  GDOP  is  oinlaua 


EXPERIMENT 


ESTIMATE 

PARAMETER 


CONSIDER 

PARAMETER 


0.2«r<0.8  1.0«j<2.2  0.4<cr<1.8 


10<o<310  10<cr<480  610<o<860 


20«rc50  20«j<100  690<a<820 


o»i<a<o.  2  o.i«j<o:2  o.i«t<o. 


N  Latitude  Coaponent 
E  Longitude  Coaponent 
U  Vertical  Coaponent 
R  Radial-track  of  satellite 
B  Clock  Bias 
D  Clock-  Drift 
A  Clock  aging 

C  9  units  in  asters) 
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ESTIMATE 

PARAMETER 

P 

(1,2) 

P 

(1,3) 

P 

(1,4) 

P 

(2,3) 

P 

(2,4) 

P 

(3,4) 

N(l) 

E  (2) 

U  (3) 

T  (4) 

-.2<p<.4 

-.2<p<— .1 

7<p<-.3 

6<p<.l 

-.5<p<.4 

. 6<p<.8 

N  Cl) 

E(2) 

U  (3) 

L(4) 

0<p<. 6 

0<p<.5 

• 3<p< • 7 

.5<p<.6 

.5<p<.7 

.5<p<.6 

N  (1) 

E  (2) 

U(3) 

C  (4) 

-.5<p<-.3 

.5<p<.7 

-. 7<p<-.5 

-.5<p<.l 

. 7<p<.8 

.7<p<.8 

N  (1) 

B(2) 

U  (3) 

*(4) 

-.2<p<.5 

— . 4<p<.l 

— .5<p<.l 

— .  8<p<— .2 

-.7<p<-.5 

,8<p<.9 

N  Latitude  Component 
E  Longitude  Component 
U  Vertical  Component 
T  Tropospheric  correction 
L  Long-track  of  satellite 
C  Cross-track  of  satellite 
a  Radial-track  of  satellite 


C  p  unitless) 


14 


-.4<p<-.2 


.2<p<.4 


-.6<p<-.2 


,6<p<.2 


-.9<p<—.7 


15 


•.3<p<e. 


-.3<p<— .2 


-.2<p<.2 


-.4<p<-.2 


.1<P<.6 


A  (4) 


N  Latitude  Component 
£  Longitude  Component 
U  Vertical  Component 
B  Clock  Bias 
D  Clock  Drift 
A  Clock  Aging 

(  p  unitless) 
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